Contiguous redshift parameterizations of the growth index 
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The growth rate of matter perturbations can be used to distinguish between different gravity 
theories and to distinguish between dark energy and modified gravity at cosmological scales as an 
explanation to the observed cosmic acceleration. We suggest here parameterizations of the growth 
index as functions of the redshift. The first one is given by 7(a) — 7(a) jtt^ — 7^ +7^ 



that interpolates between a low/intermediate redshift parameterization 7(a) = "/i^t^{a) = 70 + (1 — 
On I o,)la and a high redshift "f^^^i^ constant value. For example, our interpolated form 7(a) can be used 

. when including the CMB to the rest of the data while the form 7,^^, (a) can be used otherwise. It is 

' found that the parameterizations proposed achieve a fit that is Ijetter than 0.004% for the growth 

CN I rate in a ACDM model, better than 0.014% for Quintessence-Cold-Dark-Matter (QCDM) models, 

^! and better than 0.04% for the flat Dvali-Gabadadze-Porrati (DGP) model (with = 0.27) for the 

r-{ ' entire redshift range up to ^cmb • We find that the growth index parameters (70, 7a) take distinctive 

values for dark energy models and modified gravity models, e.g. (0.5655, —0.02718) for the ACDM 
model and (0.6418, 0.06261) for the flat DGP model. This provides a means for future observational 
■ data to distinguish between the models. 

PACS numbers: 95.36.-(-x;98.80.Es;04.50.-h 

o 

^ ; I. INTRODUCTION 

O ^, Cosmic acceleration can be caused by a dark energy component in the universe or a modification to the Einstein 
' field equations of General Relativity at cosmological scales. The growth rate of matter perturbations has been the 
subject of much recent interest in the literature as a way to distinguish between one possibility or the other, see for 
example P, S i i, H, S 0, i i, 0, [H E 111 

for a partial list. Indeed, distinct gravity theories may have degenerate 
^ I expansion histories but can be distinguished by their growth rate functions. 

As usual, the large scale matter density perturbation d — Spm/pm satisfies, to linear order, the differential equation 



S + 2HS-47rGeffPmS = 0, (1) 



where H is the Hubble parameter and the effect of the underlying gravity theory is introduced via the expression for 
\ Gef f ■ The distinct behavior of S for different gravity models can be seen in some of the aformentioned references such 
■ as for example [§, Equation ([T]) can be written in terms of the logarithmic growth rate / ~ d\n5/d\na as 



where primes denote d/d\na. Throughout this work we will use the numerically integrated solution to this equation 

• j-H normalized at a = (z = 00). Next, the growth function / is usually approximated using the ansatz [ii [n, m, m 

X ■ 

H : / = "™ (3) 



where 7 is the growth index parameter. Reference |14| made an approximation that applies to matter dominated 

models and proposed f{z = 0) = fiJJ^o ^'^'^ followed by a more accurate approximation f{z = 0) = fi^Q in p^.[l^. 
Reference jl7'| considered dark energy models with slowly varying equation of state, w, and found an expression for 7 
as function of ilm and w. This has been discussed further in more recent references, see for example [3, 19], and also 
expanded to models with curvature in ^ and plj . 

The approaches of expanding the growth index around some asymptotic value or early, matter dominated times 
with rim ~ 1, or those considering specific redshift ranges to approximate 7 do not cover other redshift ranges of 
interest where observational data is available and can constrain the growth parameters or break degeneracies between 
them and other cosmological parameters. 
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Some observational data is already available over the redshift range z = — 3.8 [1^, [H, [23, [H, [H, [13, [Ml [H, [H, [s^] 
and some recent papers have put some constraints on the values of a constant growth index parameter, see for example 
[3^ where using the cosmic microwave background (CMB), type la supernovae (SNIa), and X-ray cluster gas-mass 
fractions (/gas), the authors found 7 = 0.5i;to'i5 and — 0.274^Q Q^g (68.3 per cent confidence limits), for a flat 
ACDM background. Also, reference ^] considered early dark energy (EDE) models and combined data from WMAP 
five-year data release, baryon acoustic oscillations and type la supernovae luminosity distances, measurements of the 
linear growth factors, Gamma-Ray Bursts (GRBs) and Lyman-a forest, obtaining 7 = 0.622 ± 0.139 (Icr) error bar, 
which is in agreement with the ACDM and with the values obtained by other papers [25[, [2^, and as stated there, 
with slightly smaller error bars. Their error bars on 7 were similar to those forecasted for future weak lensing and 
SNIa data by [ll[3i]. 

In this analysis, we propose parameterizations of the growth that are function of the redshift, covering a wide range 
of low and intermediate redshifts, and then transitions after some redshift to an almost constant growth index at very 
high redshifts up to Zpj,,,^ at the decoupling epoch. 



II. INTERPOLATED PARAMETERIZATION OF THE GROWTH INDEX 

First, we recall here the work of references [3, '22'] where the authors proposed a redshift dependent parameterization 
of the growth index that was intended for the redshift range < z < 0.5 [1, [l^l and reads 

7(2:) = 70 + 7' 2: (4) 

where 7' = ^(z = 0). This showed already the potential of a variable growth index to distinguish between dark 
energy models and modifi ed g ravity models \4. [ 231 . However, current growth data is already available over the redshift 
range z = - 3.8 [H, [H, [ll [H, [2I, [13, [U,]!!!! [3 , well beyond the z <! approximation. 

In order to be able to consider constraints from this higher redshift data, future growth data that spans over higher 
redshift ranges, and also very high redshifts up to the CMB scale (whether to break parameter degeneracies or to 
put direct constraints), we propose a parameterization that covers such wide ranges and interpolates to an almost 
constant value of 7 at very high redshifts up to z^^^^ . Similar to the interpolation proposed in (23j (see appendix C 
there) for the equation of state of dark energy, we propose here the following parameterization for the growth index: 

7(a) = 7(a) /i(a/a„J -|- jeariy (1 - h{a/a^^J) (5) 

where the subscript ttc stands for transition to a constant (or almost constant) early growth index, jeariy The 
function h is chosen to have the following property 

Ha/a^J^O for a « = 1/(1 + z„ J 

'»(a/a«c)^l for a > = 1/(1 + z„J. (6) 
For simplicity, we adopt the interpolating function [2^ that achieves the behavior above and that is given by 

Ma:) = i[tanh(ln(x) + 1] = ^, (7) 
2 X + 1 



^(a) = 7(a) , , + 7„H« x (8) 



and propose the following form for the index parameterization 

1 1 

so that 7(a) interpolates between the asymptotic 7^^^,^, value at high redshifts {z ^ Zttc) and 

7(a),at. = 7(0) = 70 + (1 - a)ja (9) 
at lower redshifts (z ^ zttc)- Similarly, using a — 1/(1 + z), our parameterization reads 

7(^) = 7(2) l+, + learly (^0) 

and interpolates between pearly (note that this was also noted in the literature as 700) at high redshift (z > Zuc) up 
to the CMB scale and the following form at lower redshifts, i.e. for z < zuc 

l{.z)iate = l{z) = 70 + (y^) ^'^ 
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FIG. 1: GR - Dark Energy Models. TOP LEFT: We consider the QCDM models with a constant equation of state and 

plot the relative error ""j. in order to compare the fit of the proposed parameterization to that of the growth rate, /, 
that is numerically integrated from the growth ODE. For the ACDM, we find the best fit parameters 70 — 0.5655 and ja ~ 
—0.02710 when 7^*^^*^ = 6/11. The fit approximate the growth function / to better than 0.004% while the best fit constant 
'Jconft^' = 0.5509 approximates the growth to 0.6%. Using our redshift dependent parameterizations of growth index provides 
an improvement to the fit of the growth of about a factor 150. TOP RIGHT: We plot 7(2) = 7(2) — j^r h 700 — — 

for various values of the constant equation of state w showing very little dispersion of the order of 0.015 at any given redshift. 
BOTTOM LEFT:We consider the QCDM models with a variable equation of state, as well as some Early Dark Energy models 

and plot the relative error — - in order to compare the fit of the proposed parameterization to that of the growth rate, /, 
that is numerically integrated from the growth ODE. We find using our redshift dependent parameterizations of the growth 
index are able to approximate the growth to within 0.15%. BOTTOM RIGHT: We plot 7(2) = 7(2) — +7^0 — rr^ — 

for various dark dnergy models with a varying equation of state w{a) including some early dark energy models. 



It is worth clarifying that the zuc here is the redshift of transition from a varying growth index parameter 7(2) to 
an almost constant one, i.e. jeariy (or 700 )■ This is not necessarily the same ztrans that characterizes the transition 
from a decelerating cosmic expansion to an accelerating one. 

We show in the next sections that the proposed parameterizations fit very well the growth function that is numer- 
ically integrated from the differential equation ([2]) for a given theory. The fit is better then 0.004% for the ACDM 
model for the entire range of redshift from to the z^^^^^ — 1089 and better than 0.04% for the flat DGP model with 
fij^ — 0.27. We discuss application to these and other models in the next sections. 



III. DARK ENERGY MODELS 



For the spatially flat dark energy models with constant equation of state w, the Friedmann equations give 

§^^-l[^+w{l-nm)]. (12) 
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and energy conservation equation read 

Now, substituting Eqs. p2|) and (fTS]) into Eq. one gets 



3w;fi,„(l - f2„i) 



(film 



Next, using / = into Eq. HH) yields 

d'y 



SwQrni^ — ^m) In 17,; 



3wn^{j - 1/2) + f);; - ^nl-'T + 3w-f-^w + ^^o. 



An expansion of Eq. (fT5|) around il„j = 1 (early times), to the first order of (1 — flm), one gets [13, [13 



7 = 



3(1 -w) 



3 (l-w)(l-3w/2) 
5-6w ' 125(l-6w/5)2(l-12zi;/5) 

Eq. p6p gives the asymptotic expression for very high redshifts 



(1 - a™)- 



3(1 



w) 



5 — 6w 



(13) 



(14) 



(15) 



(16) 



(17) 



which reduces to the wch-known 7^ 
the form 



= for the ACDM model so that our parameterization in this case takes 



7(z) = 7(z) 



l+z 



1 + 



1 + 



(18) 



l+z 



with 7(2) given by our low/intermediate redshift parameterization (jlip . 

We show in Figure [1] how well various parameterizations of the growth fit the growth rate function / that is 

integrated numerically from the differential equation (|14p by plotting the relative error — -. For ACDM, we plot 
the figures up to a redshift of 5 but we performed best fits of the parameters up to the Zp^g ~ 1089. We find the 
best fit parameters 70 = 0.5655 and 7a — —0.02718 when 7^*^^^^ = 6/11. The fit approximate the growth function / 
to better than 0.004% while the best fit constant Jcon?t^ — 0.5509 approximates the growth rate function to about 
0.6%. Thus using our redshift dependent parameterizations for the growth index gives an improvement to the fit of 
the growth rate function of about a factor 150. We also plot 'y(z ) for various dark energy equations of state, including 
some early dark energy models (see for example [H, [sj [H, \^ for a discussion of the latter and our appendix for a 
parameterization). In Table |T] we list the best fit parameter values of 70 and 7a for the various models used. We find 
that these best fit values for 70 and 7a do not change for a wide range of zttc from for example 0.5 to several. 



IV. DGP MODEL 



For the spatially flat DGP [2J| model, the effective gravitational constant is given by 



and Friedmann equations yield 



The conservation of energy gives 



G Hi + ni) 



H 30, 



1 -|- 



(19) 



(20) 



(21) 



5 




FIG. 2: DGP models. LEFT: We consider the DGP model and plot the relative error ~^ in order to compare the fit of 
the proposed parameterization to that of the growth rate fnum that is numerically integrated from the growth ODE. We find 
the best fit parameters 70 — 0.6418 and ja = 0.06261 for Q'^ — 0.27 when "f^^ — 11/16. The fit approximates the growth 

function / to better than 0.04% while the best fit constant 'jconst = 0.6795 approximates the growth to 1.95%. So using our 
redshift dependent parameterization of the growth index provides an improvement to the fit of about a factor 50 for the DGP 
model. RIGHT: We plot 7(2:) — 7(2) — h 7^0 — — for various values of Jl^ showing very little dispersion of the 

order 0.01 or less at any redshift. 



where the matter energy density is given by 



" [(1 - + v/f^o^(l + z)3 + (1 - f}0,)V4]2 ■ 



Now, using equations HH), ^ and ^ into Eq. one gets 
Next, using / = fi^n into Eq. (1^ . one gets 

_ 3r!„,(i-a„)ina„ dj _ 3(1 ~ n,nh 07 . 2 ' _ r!,\rT(i + 2r>^) ^ 

1 + n„, dn„, 1 + a„, i + n„, i + n^ ' ^ ' 

Again, expanding Eq. ([24]) around Q,n = 1, to the first order in (1 — il„i), one gets Q 

So the asymptotic value for the DGP model is ^^^^ = Using the proposed parameterizations ((TO)l and (fTTjl for 
the DGP model, we show in Figure [2] how well the parameterizations fit the growth rate function / that is integrated 

numerically from the differential equation For that, we plot the relative error ^"^ . We performed the fit 

for a redshift up to the z^^^g ~ 1089. We find the best fit parameters 70 = 0.6418 and 7a — 0.06261 for fij^ = 0.27 
when 7^^^^ = 11/16. The fit approximates the growth function / to better than 0.04% while the best fit constant 
Iconst — 0.6795 approximates the growth to 1.95%. Thus, using our redshift dependent parameterizations of the 
growth index provides an improvement to the fit of the growth rate function of about a factor 50 for the DGP model. 
We also plot 7(2:) using our parameterization for various values of U,m and find that the difference of the order 0.01 or 
less at any redshift. In Table U we list the best fit parameter values of 70 and 7^ for the various models used. Again, 
we find that these best fit values for 70 and 7a do not change for a wide range of zuc fi'om for example 0.5 to several. 



V. CONCLUSION 

Data on growth rate of large scale structure covers already a wide range of redshift [1^, [H, |23, [H, [2^, [s^, HH, 
[3^ . [ssl . [s^ l and is likely to cover even a wider range for incoming and future data. We proposed parameterizations 
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Parameters for various QCDM models. 


(wo,Wa) 


70 


7a 


(-0.8,0) 


0.5690 


-0.02131 


(-0.9,0) 


0.5683 


-0.022525 


(-0.95,0) 


0.5676 


-0.02699 


(-1,0) 


0.5655 


-0.02718 


(—1.05, 0) 


0.5635 


—0.02735 


(-1.1,0) 


0.5617 


-0.02749 


(-1.2,0) 


0.5583 


-0.02771 


(-1,0.11) 


0.5641 


-0.02464 


(-0.8,-0.3) 


0.5720 


-0.03074 


(-1.2,0.8) 


0.5409 


-0.01417 


Parameters for some EDE models. 


(wo,C) 


7o 


7a 


(-0.972, 1.858) 


0.5498 


-0.02915 


(-0.95,2.5) 


0.5165 


-0.05578 


Parameters for various DGP models. 




70 


7a 


0.22 


0.6314 


0.07324 


0.27 


0.6418 


0.06261 


0.32 


0.6504 


0.05279 



TABLE L We list the parameter values for in our interpolation parameterization for various QCDM, EDE, and DGP models. 
These values were found by fitting our parameterization to the numerically integrated solution of ODE for the growth function, 
/ (e.g. we use for 7(2), Eqs.(18) with(9) for dark energy models, and Eqs. (25) with (9) for DGP models). We see that the 
QCDM and EDE models have a negative values for the parameter 7a, while the DGP models have a positive value for 7^, 
thus providing parameter that observational data can constrain to distinguish between the two gravity theories, additionally 
70 takes on distinct values for each theory. 

of the growth index that cover such wide redshift ranges and also interpolates to the highest redshifts including the 
CMB scale. The parameterizations are found to fit the growth function to better than 0.004% over the entire redshift 
range for the ACDM model, to better than 0.014% for various QCDM models, and to better than 0.04% for the DGP 
model (with = 0.27). Such parameterizations should be useful for ongoing and future high precision missions. 
We find that the best fit values for the growth index parameters take distinctive values for dark energy models versus 
modified gravity models: (70,7a) = (0.5655,-0.02718) for the ACDM model and (70,7a) = (0.6418,0.06261) for the 
flat DGP model. Most notable of the above values is the fact that is of a different sign for the two model. This 
distinction hold even when looking at more complex dark energy models. This provides a way for observational data 
to distinguish between dark energy models and modified gravity models as cause of cosmic acceleration. 
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APPENDIX 

For early dark energy we use a Mocker model first introduced by ^35i] . The dark energy equation of state for these 
models is given by: 

(26) 

In these models the dark energy component behaves like nonreletavistic matter at high redshifts, having an equation 
of state w = 0, but assymptotes to a cosmological constant with w = —1. We use parameter values for wq and C 
given by [sj which are said to fit CMB and SN la constraints very well. See the aformentioned references as well as 



w{a) 



l + Wo 
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[40| for a more in depth description of these models. 



[1] A. Lue, R. Scoccimarro, G. Starkman, Phys. Rev. D 69, 124015 (2004) 

[2] V. Acquaviva, A. Hajian, D.N. Spergel and S. Das, Phys. Rev. D 78, 043514 (2008). 

[3] Y.G. Gong, Phys. Rev. D 78, 123010 (2008). 

[4] D. Polarski and R. Gannouji, Phys. Lett. B 660, 439 (2008). 

[5] K. Koyama, R. Maartens, JCAP 601, 16 (2006) 

[6] T. Koivisto and D.F. Mota, Phys. Rev. D 73, 083502 (2006). 

[7] S. Daniel, R. Caldwell, A. Cooray and A. Melchiorri, Phys. Rev. D 77, 103513 (2008). 

[8] L. Knox, Y.-S. Song and J.A. Tyson, Phys. Rev. D 74, 023512 (2006). 

[9] M. Ishak, A. Upadhye and D.N. Spergel, Phys. Rev. D 74, 043513 (2006). 
[10] E.V. Linder, Phys. Rev. D 72, 043529 (2005). 
[11] I. Laszlo and R. Bean, Phys. Rev. D 77, 024048 (2008). 
[12] B. Jain, P. Zhang, Phys. Rev. D78, 063503, (2008). 
[13] W. Hu, I. Sawicki, Phys. Rev. D76, 104043, (2007). 

[14] P.J.E. Peebles, The Large-Scale Structure of the Universe (Princeton University Press, Princeton, New Jersey 1980). 
[15] J.N. Fry, Phys. Lett. B 158, 211 (1985). 

[16] A.P. Lightman and P.L. Schechter, Astrophys. J. 74, 831 (1990). 
[17] L. Wang and P.J. Steinhardt, Astrophys. J. 508, 483 (1998). 

[18] We find a small discrepancy in the denominator of the second term of equation (|16|l compared to the widely cited expression 

for the growth index, i.e. equation (Bll) from pioneering work by [l7j . 
[19] E.V. Linder and R.N. Cahn, Astropart. Phys. 28, 481 (2007). 
[20] Y. Gong, M. Ishak, A. Wang, arXiv:0903.0001vl [astro-ph.CO] (2009). 
[21] M. J. Mortonson, W. Hu and D. Huterer, Phys. Rev. D 79, 023004 (2009). 
[22] R. Gannouji and D. Polarski, J. Cosmol. Astropart. Phys. 05 (2008) 018. 
[23] E. Komatsu et ai, Astrophys. J. Suppl. 180, 330 (2009). 
[24] G. Dvah, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000). 
[25] C. Di Porto and L. Amendola, Phys. Rev. D 77, 083508 (2008). 
[26] S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 77, 023504 (2008). 
[27] L. Guzzo et al., Nature 451, 541 (2008). 

[28] M. CoUess et ai, Mont. Not. R. Astron. Soc. 328, 1039 (2001). 
[29] M. Tegmark el al, Phys. Rev. D 74, 123507 (2006). 
[30] N.P. Ross et ai, Mont. Not. R. Astron. Soc. 381, 573 (2007). 
[31] J. da Angela et ai, Mont. Not. R. Astron. Soc. 383, 565 (2008). 
[32] P. McDonald et al, Astrophys. J. 635, 761 (2005). 

[33] M. Viel, M.G. Haehnelt and V. Springel, Mont. Not. R. Astron. Soc. 354, 684 (2004). 
[34] M. Viel, M.G. Haehneh and V. Springel, Mont. Not. R. Astron. Soc. 365, 231 (2006). 
[35] E. Linder, Phys. Rev. D 73, 063010 (2006). 

[36] D. Rapetti, S. W. Allen, A. Mantz, H. Ebehng. lMXiv:08lOy59V l [astro-ph] (2008). 
[37] J.Q. Xia, M. Viel, JCAP 0904 002 (2009). 

[38] L. HoUenstein, D. Sapone, R. Crittenden, B. M. Schaefer, JCAP 0904:012 (2009). 
[39] D. Huterer and E. V. Linder, Phys. Rev. D 75, 023519 (2007); 
[40] E. Linder, Astropart. Phys. 26, 16 (2006). 



